Ehrenfest-time Dependence of Counting Statistics for Chaotic Ballistic Systems 
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Transport properties of open chaotic ballistic systems and tiieir statistics can be expressed in 
terms of the scattering matrix connecting incoming and outgoing wavefunctions. Here we calculate 
the dependence of correlation functions of arbitrarily many pairs of scattering matrices at different 
energies on the Ehrenfest time using trajectory based semiclassical methods. This enables us to 
verify the prediction from effective random matrix theory that one part of the correlation function 
obtains an exponential damping depending on the Ehrenfest time, while also allowing us to obtain 
the additional contribution that arises from bands of always correlated trajectories. The resulting 
Ehrenfest-time dependence, responsible e.g. for secondary gaps in the density of states of Andreev 
billiards, can also be seen to have strong effects on other transport quantities such as the distribution 
of delay times. 

PACS numbers: 03.65.Sq, 05.45.Mt 



I. INTRODUCTION 

After the conjecture by Bohigas, Gianonni and Schmit 
in 1984-'^, that chaotic systems are well described by ran- 
dom matrix theory (RMT)^, research started to demon- 
strate this connection on dynamical grounds by means 
of semiclassical methods based on analyzing energy av- 
eraged products of expressions similar to the Gutzwiller 
trace formula'^ for the density of states, that are asymp- 
totically exact in the limit ?i — >■ 0. For open systems 
we are particularly interested in the scattering matrix 
S{E), which is an TV X A^ matrix if the scattering leads 
carry A^ states or channels in total. Its elements can, like 
the Gutzwiller trace formula, be expressed'* in terms of 
sums over the classical trajectories containing the stabil- 
ity factors of the orbits A^ and rapidly oscillating phases 
depending on the classical actions S^ of the considered 
trajectories 7 divided by h 



^O'l 




A.^e^^/^)s, 
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with Tfj = 2'KhA with the mean level spacing of the quan- 
tum system A. Here the sum is over the scattering trajec- 
tories that connect the two channels i and o. For systems 
with two (or more) leads the scattering matrix breaks up 
into reflecting and transmitting subblocks, so we might 
restrict our attention to trajectories starting and ending 
in certain leads. 

In the context of spectral statistics, i.e. for the two 
point correlation function of the density of states con- 
taining a double sum over periodic orbits, this dynamical 
understanding of the conjecture^ was - as for other quan- 
tities - achieved in several steps. Starting with the pair- 
ing of identical (or time reversed orbits in the presence 
of time reversal symmetry) the so called diagonal con- 
tribution was evaluated in Ref. 5 using a sum rule from 
Ref. 6. Nondiagonal contributions consisting of pairs of 
long orbits differing essentially only in the place where 
one of the orbits possesses a self crossing and the other 



avoids this crossing were analyzed in Ref. 7. This was 
extended^ and formalized for orbits differing at several 
places, so called encounters. 

In the context of transport, i.e. for example for the 
two-point correlator of scattering matrix elements, which 
if restricted to the transmission subblocks is via the 
Landauer-Biittiker formalism^ proportional to the con- 
ductance, the diagonal contribution was calculated in 
Ref. 10. An orbit pair differing only in one crossing was 
analyzed in Ref. 1 1 and this was again extended to orbits 
differing at several places^^. These results and those for 
closed systems agreed with results from RMT, but be- 
sides this dynamical understanding of the RMT results, 
these semiclassical calculations proved very successful in 
determining the effect of a finite Ehrenfest time te on 
transport quantities, starting with the pioneering work 
of Ref. 13. The Ehrenfest time^'' separates times when 
the time evolution of a particle follows essentially the 
classical dynamics from times when it is dominated by 
wave interference. Its value is obtained as the time when 
two points inside a wave packet initially of quantum size 
h/pp with the Fermi momentum p-p evolve to points with 
a distance L of the linear system size. We thus get due 
to the exponential separation of neighboring trajectories 
in the chaotic case 



TE 



1 pfL 



(2) 



with the Lyapunov exponent A. 

Before these semiclassical calculations of the 
Ehrenfest-time dependence, there already existed 
theories to describe the effect of a finite Ehrenfest 
time on the correlators of scattering matrix elements: 
Aleiner and Larkin obtained*^ for the correlator of 
two transmission matrices, i.e. the conductance, an 
exponential suppression with increasing Ehrenfest time 
in agreement with semiclassics. This work was however 
unsatisfactory in one main aspect: a small amount of 
impurity scattering was introduced by hand to imitate 
the effects of diffraction in a ballistic system. 



Another phenomenological theory to describe the effect 
of a finite Ehrenfest time is effective RMT^^. It sphts 
the phase space and thereby also the underlying scat- 
tering matrix of the considered system into a classical 
and a quantum part, where the first one is determined 
by all trajectories shorter than t^ and the second one 
by all trajectories longer than te, as well as introduc- 
ing an artificial phase dependent on the Ehrenfest time. 
The predictions of this theory are only partially correct: 
weak localization is predicted to be independent of the 
Ehrenfest time, while the previously mentioned theories 
and also numerical simulations^^'^^ predict it to decay 
with the Ehrenfest time. In contrast to the quantum 
correction of weak localization, effective RMT gave good 
predictions for effects at leading order in N such as shot 
noise^^"^^ or the gap in the density of states of a chaotic 
Andreev billiard^^'^"*. 

Staying only at the leading order in inverse channel 
number we will consider the correlation function of 2n 
scattering matrices at alternating energies defined as 



Ci<^,ri,T) = — Tr 



^•^-^X-^) 



(3) 



where for simplicity the energy e is measured with re- 
spect to the (Fermi) energy E and in units of the so 
called Thouless energy Et — ?i/2td with the dwell time 
td measuring the typical time a particle stays inside the 
system. The latter is related to the Heisenberg time Th 
via the relation Th — Ntj:,. The Ehrenfest-time depen- 
dence is incorporated in r = te/td- The explicit form 



Ci{e,T,n) 
C2{e,T,n) 



Ci(e,T, n) + C2(e,T, n), 
C(e,n)e-^(i-'"^\ 
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(4) 
(5) 

(6) 



with the RMT (i.e. t — 0) part of this correlation func- 
tion denoted by C{e,n). The term in (5) derives from 
effective RMT^^'^^. Although this theory describes cer- 
tain phenomena quite well, e.g. the dependence of the 
Andreev gap on the Ehrenfest time^'*, a dynamical justi- 
fication of this result is still lacking. So far Ref. 25 cal- 
culated C(e, T, n) for n = 1, 2, 3 while Refs. 18,26 showed 
the separation into two terms in (4) to be a consequence 
of the preservation under time evolution of a phase-space 
volume of the system. Moreover they also calculated 
the explicit form we give in (6) for the second term and 
that the first term in (5) is proportional to the factor 

g-r(l-ine) 

Because of (1) the correlation function can be writ- 
ten semiclassically in terms of 2n scattering trajectories 
connecting channels along a closed cycle like in Fig. la. 
This leads to trajectory sets with encounters as in Fig. 
lb,c which can then be moved into the leads to create the 
remaining diagrams in Fig. 1. Including the correct pref- 
actors and the energy dependence, the correlation func- 




FIG. 1: (Color online) The trajectory sets with encounters 
that contribute to the 3rd correlation function C{e, 3) 



tion becomes semiclassically 
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T~f are the times trajectories 7 spend inside the system, 
and we identify the channels in+i = h- Note that (7) and 
this identification imply that the trajectories and their 
partners (traversed in reversed direction) considered in 
C{e,T,n) form a closed cycle. 

In this paper we want to show how (4,5,6) can be ob- 
tained using the trajectory based methods developed in 
Refs. 7,8,11,12. In Sec. II we consider the first term in 
(4): we show that the prefactor C(e, n) of the exponential 
is indeed given by the RMT expression obtained in Ref. 
27 and that this is multiplied by the exponential given 
in (5). The underlying diagrams considered here are the 
same as the ones occurring also in the semiclassical calcu- 
lation of the RMT contribution. In Sec. Ill we consider 
the second term in (4) and show how this contribution 
arises from trajectories that are always correlated. Fur- 
thermore we show in Sec. IV that there exist no mixed 
terms between the first and the second term in (4), that 
could result - expressed in terms of the considered di- 
agrams - from correlations between trajectories always 
correlated with each other on the one side and trajecto- 
ries only correlated with each other during encounters on 
the other side. 



II. INFLUENCE OF THE EHRENFEST TIME 
ON TRAJECTORIES WITH ENCOUNTERS 

The main idea in this section is to split our diagrams 
in a different way compared to the semiclassical anal- 
ysis without Ehrenfest time (referred to as the RMT- 
treatment) and the analysis of the Ehrenfest-time depen- 
dence of the cases n = 1,2,3 in Ref. 25: in the semi- 
classical calculation one considers an arbitrary number 



FIG. 2: (Color online) A 3-encounter as it can be approxi- 
mated in the RMT-treatment (c.f. Fig. Ic). The encounter 
stretches are marked by a box (shown red). 





FIG. 4: (Color online) A diagram with two 2-encounters as 
we treat it with Ehrenfest time. The encounter stretches of 
the two 2-encounters are marked by boxes (shown red and 
blue). A possible position of the Poincare surface of section 
(PSS) is marked by a black vertical line. 



FIG. 3: (Color online) A 3-encounter as previously treated 
with Ehrenfest time^^. The encounter stretches are marked 
by a box (shown red). 



Two 2-encounters 



of orbits encountering each other. It turns out in the 
RMT-treatment to be sufficient to consider only encoun- 
ters where all orbits are linearizable up to the same point, 
see for example Fig. 2. When taking into account the 
Ehrenfest-time dependence this is no longer sufficient as 
was first shown in Ref. 25; see Fig. 3 for an example of 
an additional diagram analyzed in this case. The main 
complication arising in Ref. 25 is then to treat these en- 
counters. To simplify the calculation we imagine these 
encounters being built up out of several encounters, each 
of which consist of two encounter stretches. We have dis- 
tinguished these 2-encounters by different boxes in Fig. 
4. In this way it is much easier to consider encounter dia- 
grams of arbitrary complexity with finite Ehrenfest time, 
which did not appear in the formalism used in Ref. 25. 

We first illustrate our procedure by considering three 
correlated orbits with two 2-encounters as in Fig. 4 and 
show how the result given in Ref. 25 can be obtained in 
this case and then we treat the general case of n orbits 
with (n — 1) independent or overlapping 2-encounters. 



A. Explanation of our procedure for n = 3 

In the treatment of the RMT-type contribution (5) we 
first consider the case in which all the encounters occur 
inside the system. For n — 3 we have the two semiclassi- 
cal diagrams in Fig. lb,c which include a trajectory set 
(of three original trajectories and three partners) with 
two 2-encounters in Fig. lb and a single 3-encounter in 
Fig. Ic. By shrinking the link connecting the two encoun- 
ters in Fig. lb we can see how we deform them into the 
diagram in Fig. Ic and we use this idea in our Ehrenfest- 
time treatment. 



For the calculation of contributions resulting from di- 
agrams differing in encounters we first need to review 
the notation and the important steps of the correspond- 
ing calculation in Ref. 8. An encounter of two orbits is 
characterized by the difference of the stable and unsta- 
ble coordinates Si and Ui measured in a Poincare sur- 
face of section (PSS) put inside the encounter; see Fig. 
4. In terms of these coordinates the duration of the en- 
counters is given by tend — 1/Aln (c^/ |siUi|) derived 
from the condition that the coordinates Si, Ui are only 
allowed to grow up to a classical constant c (which is 
later related to the Ehrenfest time). The weight function 
measuring the probability to find these encounters is ob- 
tained by integrating over all possible positions where 
the encounter stretches can be placed and dividing by 
the volume of available phase space (in the correspond- 
ing closed system) J7 and further by the durations of 
the encounters tonc.j to avoid overcounting the same set 
of correlated trajectories. The action difference AS' be- 
tween the orbits is in general given by a quadratic form 
of the coordinates Si^Ui determined by where the part- 
ner trajectories must pierce the PSS's to reconnect in the 
right way to form a closed cycle. For example for a 3- 
encounter one obtains^ AS = s[u'i + 52^*2 ~ ■5^2 where 
the prime denotes that the coordinates are measured in 
one PSS from the central trajectory. If we instead mea- 
sure the coordinates in two different sections, we obtain 
AS* — siui + S2U2 — S1U2 exp (— AAi) where the time At 
denotes the time the particle needs to travel between the 
two sections. This leads in the limit of well separated en- 
counters to AS « siUi + S2U2. From this and from Ref. 
8 we can draw the following conclusions for the form of 
the action difference in the case of an arbitrary number 
of (possibly overlapping) 2-encounters: In the case of k 
well separated 2-encounters we obtain for the action dif- 
ference AS ~ X]i=i Si^i- When these encounters overlap 
the action difference can differ from the last expression 
by terms exponentially damped with the time difference 



between the two sections. 

In our treatment, the overall contribution C*(e, r, 3) of 
the two 2-encounters (depicted in more detail in Fig. 4) 
is obtained by allowing the upper trajectory to possess 
a minimal length of the first 2-encounter and the low- 
est one a minimal length of the second 2-encounter. The 
middle trajectory, which passes through both encounters 
has a minimal length given by the maximum of the two 
encounter times as we allow the encounters to overlap. 
However we do not yet allow one encounter to be sub- 
sumed into the other so we also set the time t between 
the start of the first encounter and the end of the sec- 
ond to be longer than the maximum encounter time. To 
write down the semiclassical contribution of the diagram 
in Fig. 4 we sum over the number of possible classical 
orbits using the open sum rule^^. Converting the time 
integrals resulting from this rule to time integrals with 
respect to link durations, we obtain 



C\e,T,3) 




(8) 



X / dte('/'')^^e-*(i-'^)/^°, 

/max{tcnc,l,*onc,2} 

where the superscript refers to Fig. 4. We have summed 
over the possible channels, and ti with i = 1, . . . , 6 la- 
bel the links from the channels to the encounters. In 
(8) where d^s = dsids2 and d^u ~ duidu2, Si and 
Ui with i = 1,2 are the stable and unstable coordi- 
nate differences between the two parts of the trajecto- 
ries piercing through a PSS placed in the i-th encounter. 
As explained above, the action difference is given by 
AS* = siui + S2U2 — S1U2 exp (— AAi). By expanding 
the part of the exponential c^^/'^'l^^ containing this Ai- 
dependent part into a Taylor series one verifies easily that 
contributions from higher order terms than the leading 
(time independent) one are of higher order in 1/ {Xtjj) 
and can be neglected. This reasoning also holds for dia- 
grams with more than two 2-encounters. 

In the first line of (8) we can see that each integral 
over the links is weighted by its classical probability to 
remain inside the system for the time ti which decays 
exponentially with the average dwell time tu- We only 
want to consider trajectory sets where the whole diagram 
remains inside the system, as if any parts were to hit the 
lead and escape the diagram would be truncated at that 
point. With the energy dependence in (7) this gives the 
factors e~*»(i~"^)/'^D in (8). Inside the encounters however 
we have trajectory stretches that are so close that the 
conditional survival probability of secondary traversals 
is 1 and we need only consider the survival probability 
of one stretch. If that stretch does not escape then nei- 
ther will the other. The energy dependence still depends 
on the total time so that encounter 1 would lead to the 
factor Q-*''ac,i{i-2i(-)/TD ^ With the overlap, encounter 2 



would then have a more complicated exponential factor, 
but because the time t (between the two outer ends of 
the encounter stretches on the middle trajectory shown 
in Fig. 4) passes through both encounters their survival 
probability (of both stretches of both encounters) can be 
expressed as the survival probability of a stretch of du- 
ration t as in the last line of (8). The energy dependence 
instead also requires the extra traversal of the encounters 
as given by the exponential factor in the middle line of 

(8). 

Performing the integrals in the first line of (8) we have 



C\e,T,3)^jP^F\r), 



(9) 



(l-i6)6 

where we have moved all of the Ehrenfest-time dependent 
parts into the factor F^{t) with the superscript again 
referring to Fig. 4, 

re „(i/R)AS ie(te„c,l+tc„c,2)/rD 



d^sd^ 



^ ' ^cnc.l^cnc,2 



00 



X / dte-*(i-'')/^°. 

' max{icnc,l ,icnc,2} 



(10) 



Here we can also see the connection with the previ- 
ous Ehrenfest-time treatment of such a diagram. When 
t > ionc,i + icnc,2 the two cncountcrs separate (the inte- 
grals can then be further broken down into products) and 
this is the case in which the trajectories can be consid- 
ered to have two independent 2-encounters as in Ref. 25. 
Because we choose a different lower limit however, the 
contribution above also includes some of the diagrams 
previously treated as 3-encounters in Ref. 25. The rea- 
son for our choice becomes clear in the following steps. 
We first substitute t' = t — max{tcnc,i,ionc,2}, 



F\r) 



d^sd^u 



,(i/R)ASpie(t<,„c,i+t„„,2)/i-D 



^' ^cnc,l^cnc,2 



(11) 



and then substitute Ui — c/ai, Si — cxiUi and perforin 
the (Ti-integrals using the explicit form of the ienc,i = 
1/Aln (c^/ |siMi|) (for details of this calculation, see also 
Ref. 25). This results in 



F'^ir) = 16 /" dx^ 
Jo 
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(12) 



Now we substitute x'^ = XiC^ /h and obtain 
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FIG. 5: (Color online) One 2-encounter is located fully inside 
the other, corresponding to our treatment of a generalized 
version of a 3-encounter. The two 2-encounters are marked 
by boxes (indicated by different colors). 




FIG. 6: (Color online) The second of two 2-encounters now 
enters the lead so that only tc of it remains inside the system. 



Here we split the resulting expression into an h- 
independent integral (or more exactly trivially dependent 
on fi,), that exists due to the energy average that is always 
contained in our calculations, and an Ehrenfest-time or 
h dependent part with te = 1/Aln (^c^/h). This contains 
the Ehrenfest-time dependence that is expected from (5), 
so (13) already shows that the diagrams considered here 
yield the correct Ehrenfest-time dependence. 



2. A 3-encounter 

Now we consider the case in which one of the two 2- 
encounters lies fully inside the other one, which we will 
refer to as a generalized version of a 3-encounter, as de- 
picted in Fig. 5. 

For the Ehrenfest-time dependent part we have a sim- 
ilar contribution as in (10) with two differences: First t 
is best defined as the distance between the midpoints of 
the two different encounter stretches and so it can vary 
between 

1 



\t\ < - (max{icnc,l,icnc,2}-min{icnc,l,Wc,2}), 



\t\ l£ :^Konc,l ~ ionc,2| 



(14) 



Second the survival probability of the encounters is de- 
termined by the longest encounter stretch and is indepen- 
dent of t. The Ehrenfest-time dependent part can then 
be written as 
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2 |*cnc,l~^cnc,2 

(i^g-(max{tenc,l,tonc,2})(l-ie)/-rD_ 
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Performing the t integral and following the same steps as 
for (12,13), we find 



F'ir) 



16 / dx —— cos(xi) 

Jo ""^ A 

X cos (x' ) e~(™^''{~'"^i^^'"^2})(l-i«)/(-^'rD) 
^g-i£(lna:;+lna;2)/(ATD)g-r(l-3ie)_ /^g-j 



This part also shows an Ehrenfest-time dependence as 
expected from (5). Note that when performing the t- 
integral the result in this case is of course proportional 
to |tcnc,i — tcnc.2\ which coutaius, after the substitution 
from a; to x', two times the same terms linear in te with 
different signs that thus cancel each other. 



3. Touching the lead 

Up to now we have concentrated on encounters in- 
side the system, but apart from these diagrams we also 
need to consider diagrams where the encounters touch the 
opening as in Fig. Id-h. We will, as above, start by con- 
sidering encounters built up out of two 2-encounters and 
we focus here on how the calculation of the contribution 
is changed when encounters move into the lead compared 
to the treatment of encounters inside the system. As can 
also be found in more detail in Ref. 25 when encounters 
touch the lead one includes in the semiclassical expres- 
sions for encounters inside the system an additional time 
integral running between zero and the corresponding en- 
counter time, which characterizes the duration of the part 
of the encounter stretch that has not yet been moved into 
the lead. 

We consider two encounters with durations ionc.i and 
icnc,2 with the second encounter touching the opening as 
in Fig. Id and drawn in more detail in Fig. 6. As the 
second encounter enters the lead we now define the time 
t to be from the start of the first encounter until the lead 
and introduce the time tc which measures the part of 
the second encounter that has not yet been moved into 
the lead. We also separate the Ehrenfest-time relevant 
contribution F^{t) in this detailed calculation into two 
cases: in the first case (A); ienc,2 < ienc.i; we have F^{t) 
with the additional integral over the time tc 
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where the limits on the time integrals derive from the 
fact that the first encounter is not allowed to touch the 
lead (this would be included as a 3-encounter) and that 



the second must. Performing the time integrals this is 
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with the first and second term in the square brack- 
ets resulting from the upper and lower limit of the Re- 
integration. In the second case (B); ionc.2 > ^cnci! we 
obtain 
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where the more complicated limits derive from not al- 
lowing the second encounter to move further left than 
the first. After integrating we have 
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The last line comes from the upper limit of the sec- 
ond ic-integral and has the same Ehrenfest-time depen- 
dence as before and in line with (5). Likewise the up- 
per tc time limit for case A in (17) leads to the same 
dependence and we can conclude that the upper limits 
of the tc-integrations yield contributions similar to when 
the encounters are inside the system and with the same 
Ehrenfest-time dependence. The remaining (lower) lim- 
its of the time integrations in (17,19) give contributions 
possessing a different Ehrenfest-time dependence which 
however always yield zero in the semiclassical limit due 
to the fact that the corresponding terms contain no tonc,2 
in the exponentials containing td- Apart from the ac- 
tion difference, the only term depending on S2 , M2 is the 
l/^onc,2- The resulting expression is rapidly oscillating as 
a function of the energy* and thus canceled by the energy 
average. 

We can repeat this procedure for the remaining dia- 
grams in Fig. 1 and see that the contributions are deter- 
mined by the upper limits of the corresponding tc inte- 
grals. For the diagrams with a generalized 3-encounter 
(Fig. lg,h) this follows as for the 3-encounter inside the 
system but for Fig. Ie where the two 2-encounters enter 



different channels (and possibly different leads) there is 
an additional subtlety. The two encounters are still al- 
lowed to overlap, so that during the time t the stretch 
now connecting both channels can always be inside en- 
counters but the individual encounters are not allowed to 
connect leads at both ends. These additional possibilities 
are considered later, where if both encounters connect to 
the leads at both ends we actually have a band of cor- 
related trajectories (treated in Sec. Ill) and if only one 
does we have a mixed term (treated in Sec. IV). With 
this organization of the encounters we see that each di- 
agram has the same Ehrenfest-time dependence as when 
the encounters are inside the system, which is in line with 
(5). 



4- Intermediate summary 

The reasoning so far in this section proves the form 
of (5) for n = 3. First of all, we know that the result- 
ing contribution from the diagrams analyzed contains an 
overall factor e"'^'-"'^"^"^-'. Secondly, the remaining inte- 
grals are independent of h and thus independent of the 
Ehrenfest time. Thirdly, the diagrams we analyze are the 
same as the ones analyzed in the RMT-case in the first 
part of Ref. 27. As in the limit te — >■ we must recover 
that previous result, this implies that C(e, r, 3) in (5) is 
indeed given by the RMT-expression. 



5. Full contributions 

Before proceeding to the general case however, we first 
want to illustrate how our calculation can be used to 
obtain, apart from just the Ehrenfest-time dependence, 
the complete dependence on td and e. 

We therefore start for the two 2-encounters from Fig. 
4 from the last expression in (13) and perform first the 
t'-integral 

F'{r) = J^^°'d.x'^^cos(xl)cos(4) 

- max{— In a; J,- In 2;2}(1 — 2ie)/(ArD) 
^gmin{-lnXi,-lna;2}ie/(ATD)g-T(l-3ie) /2J^') 

where it is simpler to rewrite the result in terms of the 
maximum and minimum value of In x'^ . For calculating 
the x^-integrals we perform partial integrations (integrat- 
ing each time the cos functions) and then perform the 
resulting integrals from zero to infinity 
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16ie (1 - 2ie) 
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i U J^ 1 tt/Q 



— max4 — In a:', ,— In a^g |(1 — 2ie)/(Ar]3) 
min{- Ina:'^,- In a:2}ie/(ATD)p— r(l — 3ie) 



ie (1 - 2i£) (i_3i^^ 
rDT2(l-ie) 



(22) 



In the first line the additional terms due the partial in- 
tegration are either zero or cancel due to the energy av- 
erage. The final result in the last line of (22) can be also 
obtained by replacing max{— Inx'^, — Inoij} /A — y\ and 
minj— Ina;']^, — lna;2} /A = 2/2 and performing the inte- 
grals with respect to yi from zero to infinity. 

To evaluate the contribution from the generalized 3- 
encounter in Fig. 5 we again perform two partial integra- 
tions in (16) and obtain 
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(23) 



where we have also left out the terms from the partial 
integrations that cancel due to the energy average. 

With these results we can now show how they connect 
to the RMT-type results. For this we need to split our 
diagrams differently and first we need the result for an 
ideal 3-encounter as depicted in Fig. 2 whose contribution 
was calculated^^ to be 



F^(t) 



(1 - 3ie) (i_3i,^ 



(24) 



With the extra factors in (9) it is clear how in the limit 
te = this reduces to the RMT-type result for a 3- 
encounter as in Ref. 27. All the remaining contributions 
should be collected together as two 2-encounters, and as 
the ideal 3-encounter is included in our generalized 3- 
encounter we first subtract (24) from (23) 



F\t)-F\t) 



2tMe- 



r(l-3ie) 



(25) 



Before we add the result from our separation of two 2- 
encounters in (22) we remember that in the treatment 
we enforce that the first encounter is to the left of the 
second. The result in (25) does not have this restriction 
so we divide by 2 to ensure compatibility and then add 
the result in (22) to obtain 



F\r) = 



1 



(l-2ie) ^_^(i_3k) 



tdT2 (1-ie) 



(26) 



This then reduces to the RMT-type result for trajectories 
with two 2-encounters when te = as in Ref. 27. The 
agreement of these results with the previous Ehrenfest- 
time trcatment^^ can be seen as the result in (26) includ- 
ing both the result from two independent 2-encounters as 
well as most of the contribution of the diagram referred to 
as a 3-encounter in Ref. 25. When splitting the contribu- 
tion in a different way as in Ref. 25 this also leads to terms 
in both classes that contain different Ehrenfest-time de- 
pendencies that only cancel when summed together. 




FIG. 7: (Color online) A ladder of consecutive 2-encounters. 
The encounter stretches are marked by boxes (shown in dif- 
ferent colors). 



B. All orders 

Although up to now we have just reproduced results 
from Ref. 25, the procedure used here has the advan- 
tage that it yields a simple algorithm for determining the 
Ehrenfest-time dependence of the corresponding contri- 
butions to Ci(e, r, n) at arbitrary order. For our example 
of 71 = 3 we showed how it was possible to split the di- 
agrams into two classes that both showed the Ehrenfest- 
time dependence as expected from (5). We want now to 
show how to generalize our way of splitting considered for 
three trajectories to diagrams containing n trajectories. 



1. Ladder diagrams 

We start again with the situation in which all of the 
encounters are inside the system and by considering a 
case analogous to Fig. 4, but now involving n instead of 
three trajectories. We first take a diagram that consists 
of a ladder of (n— 1) 2-encounters so that the central n — 2 
trajectories each contain two encounter stretches while 
the two outside trajectories only contain one encounter 
stretch each. This situation is depicted in Fig. 7 and 
the encounters are thus characterized by (n — I) s, u- 
coordinates. 

In this case we obtain for the Ehrenfest-time relevant 
contribution F'^ (r) that the t-integral measuring the time 
difference between the end points of the two encounter 
stretches on the middle orbit in (10) is replaced by n — 2 
integrals over times ti with the same meaning as t; they 
measure the time difference between the end points of 
the two (consecutive) encounter stretches on the central 
trajectories containing two encounter stretches. These 
times likewise run from the maximum of the correspond- 
ing encounter times to infinity. The survival probability 
is determined by a single (artificial) stretch that runs 
through all the encounters so that the exponential term 
describing the td- and e-dependence is now given by 

(27) 
where tend are the durations of the (n — 1) individual 2- 



encounters and the middle exponential compensates for 
the fact that the middle encounters are traversed by two 
ti and that only one traversal should contribute to the 
survival probability. Setting i'^ = t^ — max{tcnc,i,ionc,i+i} 
and repeating now the steps of (12,13) we find the 
Ehrenfest-time dependent factor in this case to be 
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(28) 



again confirming the Ehrcnfcst-time dependence of (5). 



2. Single encounter 

Along with the case in which none of the encounters 
in the ladder can move completely inside another, we 
can look at the opposite extreme where all the encounter 
stretches lie inside of the encounter k with the longest 
duration tonc.fe = max^ {icnc,i} where tcnc,i are the dura- 
tions of the (n—l) individual 2-encounters with one of the 
two orbits containing the stretch of duration ienc,fc- This 
situation is like a generalization of the diagram in Fig. 5 
and we similarly now define the times ti to be measured 
between the centers of encounter i and the encounter 
k of maximum length (with i ^ k). Here the same 
Ehrenfest-time dependence 6""^'^^™^' follows by taking 
into account that each time ti has a range of variation 

of size tr-nr 



tcnc,i and that the td- and e-dependent 



exponential in this case is 



g-to„c,fe(l-ie)/rDgieEr=l'*o 
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This yields for the Ehrenfest-time dependent factor 

n— 1 n—l 

\ n ) 
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confirming again the Ehrenfest-time dependence pre- 
dicted by (5). 



3. Mixture 

Of course it is additionally possible to have a mixed 
form between these two extreme cases. This means that 
some 2-encounters only overlap like in the case of a ladder 
diagram while the others form 'single' encounters; see 




FIG. 8: (Color online) One possible example of a mixed case: 
One encounter is fully contained inside an other, the others 
form a ladder as considered before. 



Fig. 8 for a possible diagram. We then have a ladder 
of 'combined' encounters that themselves can be made 
up of one or more 2-encounters. The treatment of such 
diagrams is very similar to the treatments above, and 
the only slight complication is in defining the appropriate 
times to extract the Ehrenfest-time dependence. 

We recall that the first and last trajectories only pass 
through one 2-encounter while the central n— 2 trajecto- 
ries pass through two. Numbering the central trajectories 
from 1, ... ,n — 2, so that trajectory i has encounters i 
and i + 1 along it, we divide them into two sets: those 
whose encounter stretches lie fully inside each other or a 
connected encounter, as in the case of a single encounter 
above, that we place in the set fli. We place the re- 
maining orbits with two stretches separated as in ladder 
diagrams in the set fl2- As mentioned above, we con- 
dense the overlapping encounters into combined encoun- 
ters and record in the set A the labels of the trajectories 
that pass through the second stretch in each combined 
encounter. We also include in this set combined encoun- 
ters made of a single separated (ladder) 2-encounter. We 
then use m{i) for i e A to record the number of addi- 
tional consecutive trajectories involved in the same en- 
counter, so that m{i) = for separated 2-encounters and 
m{i) > for larger encounters corresponding to the sin- 
gle encounter case above. If the last combined encounter 
is a 2-encounter its second stretch is traversed by the last 
trajectory in the diagram which we number by n — 1 and 
include as an element of A. For example, for the dia- 
gram in Fig. 8 we would have ^i — {1}, ^2 — {2,3}, 
A = {1, 3, 4}, m(l) = 1, m(3) = and m(4) = 0. For the 
elements i € A we also label by ki the corresponding en- 
counter of maximum length among those from encounter 
i to encounter i + m{i). To be precise, the two stretches 
that stay together longest have length ioncfe; while the 
other encounter times are defined by how long the re- 
maining stretches remain close to one of the two longest. 

For the trajectories passing through two separated con- 
densed encounters i € Q2 we define the times ti to in- 
clude the whole of the leftmost and rightmost condensed 
encounters, i.e. to include the encounters kj and /ci+i, 
where i is the largest element in A that is < i. In this 
case the td- and e-dependent exponential can be written 




FIG. 9: (Color online) A general diagram containing also or- 
bits with more than two stretches. The encounters are marked 
by boxes (shown in different colors). 
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^2-* = ► 



FIG. 10: (Color online) Definition of the times ti in the case 
of more than two encounter stretches on one orbit. The en- 
counter stretches are shown thicker (blue). 
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(31) 



where A is A with its largest element removed so that the 
second term accounts for the overlap between the <i's and 
the third term for the fact that there is also no overlap 
at the start of the first such stretch. This equation incor- 
porates both (27) and (29). If we introduce the notation 
M[i, j] = max {ienc,i, iencj }, wc can then define the times 
t'^ — ti — M[fcj, fci+i] as before (28). Making the substi- 
tutions as done previously yields for the Ehrenfest-time 
dependent factor 

^'(-) - i-^) n/ d.;.cos(.;.)n / dt^ 
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(32) 



with M [i , j] = max | — In x^ , — In a;' } . As 572 and A must 
have the same number of elements, this again shows the 
predicted Ehrenfest-time dependence. 



4- General encounters 

Up to now we restricted our discussion to diagrams 
in which each trajectory is involved in one or two en- 
counters. This is however not yet the most general case 
where the only restriction is that each trajectory contains 
at least one encounter stretch, so that some trajectories 
can also contain more than two encounter stretches. Note 
that the situation where two trajectories interact (pass 
through the same 2-encounter block) more than once can- 
not occur at leading order in inverse channel number. An 



example of a diagram that is possible is depicted in Fig. 9. 
In the most general case we define the times ti slightly dif- 
ferently: first we separate the k > 2 trajectories that have 
one encounter stretch from the remaining n—k that have 
more than one. Then we number our encounters accord- 
ingly, first those along the trajectories with one encounter 
stretch with duration tcnc,i, * = 1, . . . , fc then the remain- 
ing encounters with duration tonc,i, i = fc + l,...,n— 1. 
For the n — k trajectories with two or more encounter 
stretches we now define ti, i = 1, . . . ,n — fc, to be the 
time difference between the outer edges of the outermost 
encounters along those trajectories. 

For any trajectories with more than two encounter 
stretches we will need additional time differences to fully 
fix the positions of the encounters. Because we defined 
the times ti to go through the outmost encounters, im- 
portantly the exponential factor with the survival prob- 
ability and the energy dependence does not depend on 
these additional time differences and is given by 



■ EITi'" ti{l-ie)/Tn j:"^k + l tcnc/TD ie Etl tc„o,>/rD 



(33) 



where the middle term ensures that the survival proba- 
bility only includes one copy of each encounter and the 
energy dependence involves all traversals of all the en- 
counters. 

For the remaining times we notice that, starting with 
the ladder system with 2 trajectories containing one en- 
counter stretch and n — 2 trajectories containing two 
stretches, every time we increase the number of trajec- 
tories with one encounter stretch we simultaneously in- 
crease the number with more than two. Therefore there 
are fc — 2 additional time differences needed to fix the po- 
sitions of the central encounters along trajectories with 
more than two and we define times ti for i = 1 , . . . , fc — 2 
from the left hand side of one encounter stretch to the 
right hand side of the next encounter stretch following 
on the right on those trajectories, see also Fig. 10. As 
the encounters are ordered, they are not (yet) allowed to 
be subsumed by each other or pushed past the outside 
encounters. The ranges of the times ti are then fixed by 
these restrictions. Using again M[i,j] defined after (31) 
in the following to make the notation more compact, we 
obtain for a trajectory containing m encounter stretches 
of durations tenc,i, i = 1, . . . , m, as illustrated in Fig. 10, 
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(34) 



In the second hne we substituted ? 



tj-M[j,j + l]. 



The time differences ti , which are more important for the 
Ehrenfest-time dependence, must instead just be longer 
than the maximal length of the encounter stretches lying 
on the considered trajectory. In general the numbering 
of the encounters and time differences can be more com- 
plicated than in Fig. 10 so we define l{i) to be a list 
of length m(j) of the encounters enclosed by the time 
U (including the outer encounters) and L(i) a list of 
the corresponding m{i) — 1 times t between the ends of 
those encounters. Now we can make the substitution 
t\ = ti — maXjg;(-j-) {tcnc.j}- After this substitution we 
recognize that (34) has become independent of h or the 
Ehrenfest time. Following then the steps in (12,13) we 
obtain 
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(35) 



with — In x'^-^^^ j = inaxj^m^ | — In x', } linked to the dura- 
tion of the longest encounter stretch contained within ti. 
We have also used M[i,j] defined after (32) and dropped 
the explicit i dependence of I, L and m above. Again we 
obtain the Ehrenfest-time dependence predicted by (5). 
As in the case of the ladder diagram above, we can also 
have the possibility of some encounter stretches being 
contained in larger encounter stretches and some sepa- 
rated from those larger encounters, see Fig. 11 for an 
example of a possible diagram. This just implies that 
some of the ti integrals have to be treated as was done 
in the case of the configuration shown in Fig. 5, and the 
Ehrenfest-time dependence predicted by (5) also follows 
in this case. 



5. Touching the lead 

When the encounters are allowed to enter the lead we 
again have to consider times representing how far each 




FIG. 11: (Color online) An encounter diagram containing also 
orbits with more than two stretches. In contrast to Fig. 9 en- 
counter stretches here are allowed to be contained fully inside 
others. The encounters are marked by boxes (shown in dif- 
ferent colors). 



encounter has moved into the lead (actually how much 
of the encounter remains inside the system) . As for the 
case treated in detail for n = 3 it is only the upper limit 
(namely the full encounter time) of these time integrals 
that have the necessary encounter time dependence to 
contribute in the semiclassical limit. The reasoning for 
n — 3 can then be carried over directly to the more gen- 
eral cases as the upper limits of these integrations yield 
contributions that are (up to constant factors) the same 
as the ones obtained when the encounters are inside the 
system. We thus obtain the same Ehrenfest-time depen- 
dence from encounters moved into the leads. 



C. Summary 



The separate diagrams considered in the RMT-type 
semiclassical treatment^^ can be created from the origi- 
nal collapse of trajectories onto each other and by sliding 
the individual encounters together or into the leads. The 
Ehrenfest-time treatment however suggests treating all 
of these possibilities instead as part of continuous fami- 
lies. What we have shown above in this section is that, 
if we partition this family in a particular way, for any 
partition we can find a suitable set of coordinates that 
allows us to transform the semiclassical contribution so 
that we can extract the overall Ehrenfest-time depen- 
dence. Though the exact details of this transformation 
depend on the structures of the partition, the algorithmic 
routines described above all lead to the same Ehrenfest- 
time dependence. Each partition and hence family then 
has the factor e~'^*^^~"") and no other Ehrenfest-time or 
h dependence. As we know that we must recover the 
RMT-type result C(e, n) in (5) when t^ — (since we 
treat the same diagrams) with no further Ehrenfest-time 
dependence, we then obtain the full result in (5) and 
hence provide a semiclassical justification of the effective 
RMT ansatz. 
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FIG. 12: (Color online) Band of n = 3 correlated trajectories. 
The length of the orbits is marked by a box; the duration of 
the encounter tone = 1/Aln (c^/maxi \si\ maxj |uj|) is marked 
by a dotted box. 



III. TRAJECTORIES ALWAYS CORRELATED 




^iPSS 



In this section we determine the so caUed classical con- 
tribution in (6). To obtain this contribution C2(e,T, n) 
semiclassically we consider a band of n trajectories that 
are correlated (inside the same encounter) for their entire 
duration between entering and leaving the system as in 
Fig. 12. This implies that all the trajectories have the 
same length t and that the maximum of the differences 
Si,Ui between their stable and unstable coordinates lies 
below the constant c (related to the Ehrenfest time). For 
the case n = 2 this configuration was first considered in 
Ref. 22 and then extended to n = 3 in Ref. 25. For our 
calculation we follow Ref. 25 and place a PSS at a dis- 
tance ii from the left end of the trajectories while the 
remaining time on the right of the section is denoted by 
t2 = t — ti. The semiclassical contribution C2(e,T, n) can 
be written as 

C2(e,T,n) = — / dti / dt2- 

'^D Jo Jo 
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X / as 
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(36) 



where we only include one traversal of the band in the 
survival probability and the restrictions on the s and u 
integrals ensure that the band always remains together 
under the exponential divergence of the trajectories due 
to the chaotic dynamics. Performing an integral over 
ti — ^2 and the Ui integrals gives 
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n— 1 , 



where Xi = e ^^^Si/c. Using that in the semiclassical 
limit 

d.T-sinf— j ==— e(TE-i), (38) 
with the Heaviside theta function 0{x), we finally obtain 

1 _ g-T(l-me) 
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(39) 



proving the Ehrenfest-time dependence of the C2(e, r, n) 
in (6). 



FIG. 13: (Color online) An example of a band of three tra- 
jectories that possesses an encounter with another trajectory. 
The band is marked by a thicker box (red stretches), the en- 
counter of the other trajectory with the band by a dotted box 
(blue stretch) and the links by thin (black) lines. The link 
durations are denoted by ts and t4, and the durations of the 
band before and after the PSS by fi and t2. 



IV. MIXED TERMS 

Finally we want to consider possible correlations be- 
tween trajectory structures giving the RMT-type con- 
tribution and those giving the classical part, i.e. contri- 
butions from correlations between bands of trajectories 
(that are always correlated with each other) and trajec- 
tories that are only correlated with each other during 
encounters. In particular we want to show that diagrams 
that have a correlated band that has any encounter with 
other trajectory structures (with encounters) give no con- 
tribution in the semiclassical limit. This, once general- 
ized, then excludes the existence of mixed terms in (4) so 
that (4) is complete. First we consider the case in which 
n — 1 of the trajectories form a correlated band with the 
remaining trajectory meeting the band in an encounter 
inside the system as depicted in Fig. 13. This contribu- 
tion C^^{e, r, n) to the correlation function C(e, r, n) can 
be written by treating the correlated band as before and 
introducing the times is and ^4 to represent the durations 
of the parts of the trajectory that encounter the band on 
the left and on the right of the PSS. It reads 
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where tenc is the time during which the remaining trajec- 
tory encounters the band. In this case no overcounting 



12 



factor 1/ienc occurs as the limits of the s, u-integrals in 
the third hne are chosen such that the time integral ac- 
counting for all possible positions of the single encounter 
stretch with respect to the band cancels this factor. 
Choosing different limits for the s,u-intcgrals wc could 
replace the current overcounting factor of 1/ (ti + 12) in 
(40) again by the usual 1/tcnc as in Sec. II. 

The integrals over the (n — 2) s,u-coordinates within 
the band are performed in the same manner as in the last 
section in (38) and lead to the same Heavisidc function 
0(rE — t). For the integrals in the third line in (40) over 
the differences s, u between the coordinates of a band 
trajectory and the trajectory encountering it we obtain 
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s/c. This Heaviside function is opposite 
to the one from the band, so that the total contribution 
vanishes due to the opposing restrictions of the theta 
functions. Note that when evaluating the Ehrenfest-time 
dependence of the RMT contributions in Sec. II we chose 
a different way for performing the integrals. In Sec. II we 
performed restricted time integrals first and unrestricted 
x' integrals afterwards. The different order of calculating 
these integrals leads to the theta functions which make it 
clear that mixed terms like in Fig. 13 vanish; this is why 
we use this ordering here. If we move more trajectories 
from the band (composed of at least two trajectories) to 
the trajectory structure with encounters we still obtain 
these opposing Heaviside functions and hence no contri- 
bution. 

A similar reasoning can be applied if the encounter of 
a trajectory (or part of a trajectory structure) with a 
band does not happen inside the system but enters the 
lead at the beginning or the end. In this case we obtain 
an additional time integral with respect to the time of 
the encounter that remains inside the system but, as the 
s, u-intcgrals still yield the same Heaviside functions, this 
contribution also vanishes. Note that if we move both 
ends of the encounter into the leads then the encountering 
trajectory can be considered as part of the band and 
treated as above or in Sec. HI. 

The reasoning in this section applies to an arbitrary 
number of bands of correlated trajectories connected 
by trajectories that are only correlated in encounters. 
Therefore all such mixed terms vanish. 



V. IMPLICATIONS FOR TRANSPORT AND 
SCATTERING 

A. Moments of transmission 

Up to now we have concentrated on energy dependent 
correlation functions involving the whole scattering ma- 
trix. Because we use the same semiclassical diagrams our 
result can be applied directly to dc-transport properties 
of chaotic systems such as the moments of the trans- 
mission or reflection eigenvalues. Assuming the system 
has two scattering leads and taking just the transmission 
subblock t of the scattering matrix connecting the A^i 
channels in lead 1 to the N2 channels in lead 2 (without 
an energy difference) the moments of the transmission 
eigenvalues can be written as 
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where A^ = A^i -|- A^2 is the total number of channels. 
Within the Landauer-Biittiker^ approach to quantum 
transport the moments carry information about statis- 
tical properties of the transmission in the phase co- 
herent regime, the counting statistics^^. For example 
the first moment characterizes the average conductance 
G ex M{t, 1) and the second one the power of the shot 
noise P ex M(t, 1) — M(t, 2). Using the semiclassical 
results from this paper we can simply write 



M(T,n) = M(n)e"^-h 
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where we have included the probability Ni/N of starting 
in lead 1 in the moments. Equation (43) can be gener- 
alized to ac-transport considered in Ref. 29 by including 
in the latter equation the e-dependent factors given in 
(5,6). The result in (43) again splits into two parts with 
the first involving the semiclassical moments M{n) calcu- 
lated in Ref. 30, which in turn lead to the random matrix 
probability distribution"^^'^^. 

The second term in (43) leads to the classical Bernoulli 
distribution where the transmission amplitude T is 1 
with probability N2/N and otherwise (i.e. Ni/N). The 
Ehrenfest time then provides a smooth interpolation be- 
tween these two distributions giving a weight e""^ to the 
RMT one and the remaining weight (1 — e^"^) to the 
classical one. A similar formula and result follows for 
the moments and probability distribution of the reflec- 
tion eigenvalues whose zero Ehrenfest time contributions 
can be simply derived from the treatment in Ref. 30. 



B. Moments of delay times 

Taking the full correlation functions C(e, r, n) it is pos- 
sible to obtain not only the Ehrenfest-time dependence of 
the density of states of chaotic Andreev systems, covered 
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in detail in Ref. 27, but also the moments and distri- 
bution of the Wigner delay times. We start with the 
Wigner-Smith niatrix^"^ 



Q-isHE/^. 



(44) 



which can be shown to be Hermitian by using the uni- 
tarity of the scattering matrix. Because of this unitarity 
Q can also be written as 



12. A 
Td^ 



^t 



eh 



2td7 



S 



eh \ 
2^) 



£ = 



(45) 



where the scattering matrix energy differences are mea- 
sured with respect to the energy E. The delay times are 
simply the eigenvalues of Q so their moments are 



m(T,n) = lTr[0]" 



Using the relation 



1 d" 



;r[/W-/(o)r _ =[/'(0)] 



n! de^' 

the moments of the delay times are"^"* 
d" 



m{T,n) 



'D 



i"n!Af de" 
eh 
2^ 



-Tr 



(47) 



(48) 



S^ 



S 



eh 
2^ 



e=0 



Expanding (48), the moments can be expressed as 
T?i d" 



m{T, n) 



i"n! de" ^-^ 
fc=i 



V-(_^)«-fe 



C(e,T,fc) 



e=0 



(49) 

in terms of the correlation functions calculated before. 
As this is additive we can look at the two parts of the 
Ehrenfest-timc dependent results in (4) separately. 

For the first part of the Ci (e, r, k) however, it is possi- 
ble to put our Ehrenfest dependence into the framework 
of Ref. 34 where the moments and the probability dis- 
tribution p(tw) of the Wigner delay times rw were cal- 
culated (without the Ehrenfest time). We can actually 
obtain the result in a simple way. Using (47) again we 
can see that since 



Te/= —-7- 

1 de 



5t 



2td / V 2td / 


7 
£-0 




(50) 



we have 

1 

iV 



Tr[0-TE/]" = 



d" 



i"n! de" ^-^ 
fc=i 



(-1) 



n — k 



C{e,r, k)e 



—ike 



(51) 



e=0 



Plugging in our result for Ci{e,T,k) from (5) the en- 
ergy dependent exponentials cancel so, apart from the 



damping factor e""^, we just have the moments without 
any Ehrenfcst-time dependence, leading'^'* to the RMT 
result*^^. Of course on the left hand side of (51) we have 
a simple translation by the Ehrenfest time, meaning that 
the translated probability distribution is the same as the 
RMT-one (damped). For the full Ehrenfest-time depen- 
dent distribution we simply translate back again and have 



^ / ^ v/(r+-rw)(Tw-T_) 

r± = (3±V8)TD+rE. 



(52) 



For the second 'classical' contribution in (4) we first 
take the simplest part of the contribution 



^(1) 



C^^^(e,T,fc) = 



1 



1-ifce' 



(4b) and substitute into (49) obtaining 



,(1) 



(^,^) = t5^(-i) 



fe=l 



fc" 



rSn! 



(53) 



(54) 



These moments clearly come from an exponential distri- 
bution, so that in the limit r — > oo we recover, for the 
probability distribution p2 (''^v), the classical exponen- 
tial decay of trajectories inside the system 



p«(rw) = — e'^"/^°, rw>0. 



TD 



(55) 



For the remaining contribution of the second part 
P2 (''w), we have the damping factor c^'^, and the en- 
ergy dependent phase again just leads to a shift in the 
exponential distribution so this contribution starts at te, 
thus yielding 



P2 (nv) 






TW > TE. 



(56) 



The minus sign however means we truncate the previous 
exponential at te, so the total second contribution to the 
probability distribution is 



P2(tw) 



1 



TD 



-Tw/td 



< Tw < Te, 



(57) 



and elsewhere. Since this contribution to the delay 
time probability distribution is considered to be made 
up of the two distributions pl (tw) and p\ (tw) and 
as the shifted and damped one p\ (tw) has the same 
mean (td + te) and weight as the shifted and damped 
RMT distribution /Oi(tw) but a minus sign, it is clear 
that the average time delay stays at td (i.e. from the 
untruncated exponential distribution) and is unaffected 
by the Ehrenfest time. In fact this is an example of the 
general relation^^, derived from the unitarity of the scat- 
tering matrix, in which the mean time delay depends just 
on the average spacing of the resonant levels of the scat- 
tering system and the number of scattering channels, i.e. 
it should have no Ehrenfest time or other dependence. 
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FIG. 14: (Color online) The probability density of the Wigner 
delay times for r = 1. 



This shape of the distribution however is significantly 
affected by the Ehrenfest time, and as an example we 
plot the complete probability density of the Wigner de- 
lay times for r = 1 in Fig. 14. There we can see the ex- 
ponential decay truncated at te = td before a hard gap 
separates it from the damped RMT distribution which is 
shifted to the right by te- As the RMT distribution in 
(52) starts at r_ which is (3 — V8)tj:, above the Ehren- 
fest time, the hard gap is a constant 0.172rD wide. The 
RMT distribution is continuous (unlike the truncation of 
the exponential at rg) and peaks very quickly (0.056td) 
after the gap at around 1.72e~'^/rD. Related to the semi- 
circle distribution, it likewise has an upper bound (t+), 
which interestingly can be connected"^^ to the excitation 
gap in the density of states of Andreev billiards. By ex- 
panding, for small energies, the determinantal equation 
that governs this excitation gap in terms of the Wigner- 
Smith matrix Q, Ref. 37 showed that the gap's width 
is approximately the inverse of the maximum delay time. 
There the maximum was increased by disorder scattering, 
leading to a decrease in the Andreev gap, while here the 
maximum delay time is increased by the Ehrenfest time 
as the RMT distribution is correspondingly shifted right. 
A simultaneously decreasing Andreev gap fits with the 
effective RMT^'^'^'* and semiclassical^^ treatments. How- 
ever the connection between the secondary gaps, which 
appear in the density of states of Andreev billiards^^ 
as a consequence of the Ehrenfest dependence in (4,5,6) 
shown in this paper, and the gap in the density of the 
delay times is not so clear. 



VI. CONCLUSIONS 

In this paper we have shown how to treat the effect of 
the Ehrenfest time on correlation functions of arbitrarily 
many pairs of scattering matrices. In our semiclassical 
approach we extended and combined the zero Ehrenfest 
time approach^^ (which leads to the RMT result) and the 
n — 3 Ehrenfest time approach^^ and showed how the re- 
sults of the effective RMT ansatz can be obtained. The 



different contributions are described by simple diagrams 
and following an innovative way of partitioning these dia- 
grams we implemented an algorithmic procedure that al- 
lows one to easily obtain the Ehrenfest-time dependence. 
Interestingly this always led to the same factor (which 
can be traced back to the survival probability only de- 
pending on one traversal of each encounter) so that the 
RMT-type expression is simply modified by the Ehren- 
fest time by the additional factor e""^^^"™*^). This is in 
line with the effective RMT result, but as our result is 
derived just from the underlying chaotic dynamics of the 
system we can justify for this situation the use of effec- 
tive RMT which instead conjectures the Ehrenfest-time 
dependence. 

As the semiclassical framework is based on the under- 
lying classical dynamics we can equally well move away 
from the RMT arena and obtain the 'classical' contribu- 
tion to the correlation functions. This can be seen to 
come from bands of trajectories that remain correlated 
with each other for the entire duration of their stay inside 
the system. Furthermore the fact that no mixed (between 
the RMT-type and classical-type) terms arise is simply 
due to their opposing classical restrictions. This lack of 
mixed terms as well as the classical contribution were 
previously shown to be more generally due to the preser- 
vation of volume under the dynamical evolution and the 
separation of phase-space into two essentially indepen- 
dent subsystems^^'^^. 

The separation of the correlation functions into two 
contributions, which each have a straightforward depen- 
dence on the Ehrenfest time was previously shown to be 
responsible for secondary gaps in the density of states 
of Andreev billiards^^, but has an equal effect on other 
transport quantities. For the transmission eigenvalues 
(and their moments) with no energy dependence we just 
get a straightforward interpolation between the RMT'^" 
and classical values. For the distribution of the Wigner 
delay times we further see a truncation of the classical 
(exponential) distribution and a shift to higher times of 
the RMT-type distribution. Between the two though a 
hard gap remains. 

The method described in this paper allows for the com- 
putation of the Ehrenfest-time dependence of the trace 
of arbitrarily many scattering matrix pairs but only to 
leading order in inverse channel number 1/N. The cal- 
culation was only doable because at this order the corre- 
sponding semiclassical diagrams involve no closed loops 
and have no periodic orbit encounters (surrounded pe- 
riodic orbits). When such surrounded periodic orbits 
are involved, for example for the conductance variance'^* 
or the next to leading order quantum correction to the 
transmission, reflection and the spectral form factor'^^, 
the relatively simple cancellation mechanism observed in 
this paper no longer holds. But by taking into account 
all possibilities for partial correlations within the 'fringes' 
as in Refs. 38,39 in a systematic way, it should however 
also be possible to extend our Ehrenfest-time results to 
infinite order. 
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